Emergence of geometrical optical nonlinearities in photonic crystal fiber nanowires 
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We demonstrate analytically and numerically that a subwavelength-core dielectric photonic 
nanowire embedded in a properly designed photonic crystal fiber cladding shows evidence of a previ- 
ously unknown kind of nonlinearity (the magnitude of which is strongly dependent on the waveguide 
parameters) which acts on solitons so as to considerably reduce their Raman self- frequency shift. An 
explanation of the phenomenon in terms of indirect pulse negative chirping and broadening is given 
by using the moment method. Our conclusions are supported by detailed numerical simulations. 

Photonic nanowires (PhNs), i.e. dielectric waveguides with a sub- wavelength core diameter, tight mode confinement 
and strong waveguide dispersion, have recently attracted a growing interest due to accessibility of new fabrication 
^ ■ techniques for a large variety of materials, which may lead to a number of miniaturized, high-performance photonic 
devices [1]. The small effective modal area exhibited by PhNs, which increases considerably the Kerr nonlinear 
coefficient, and the degree of controllability of the dispersion characteristics, make PhNs especially suitable for the 
investigation of extreme nonlinear phenomena such as supercontinuum generation, as many optical solitons can be 
excited by using small pump energies [2]. 

In recent theoretical work, a novel propagation equation that accurately describes the nonlinear evolution of light 
pulses in PhNs was introduced Q, see also Eq. (j2j) of the present Letter. The fundamental feature found in [J is 
that, thanks to the fact that the correct equation takes into full account of variations in the linear mode profiles of 
the waveguide with wavelength, new nonlinear effects arise in PhNs, unknown in previous formulations based on the 
generalized nonlinear Schrodinger equation (GNLSE) [4], which all assume fixed transverse field profiles. Most of the 
additional terms described in Ref. [3| have been found to have an extremely small magnitude, so that they can be 
<y5 ■ safely neglected for large core fibers. However, the longitudinal component of the electric field of the fundamental 
mode of PhNs becomes progressively more important when decreasing the core diameter or when increasing the 
refractive index contrast between core and cladding [3[. This vigorously breaks the rotational symmetry of the mode, 
whose transverse profile becomes very sensitive to frequency, thus making the new geometrical nonlinearity - as we 
shall call it in the following - extremely important. In fact, for small enough core sizes, such new nonlinear terms may 
even enter into strong competition with the Raman effect term for some range of frequencies. 

Only the circular geometry for PhNs has been considered in Ref. [3]. Circular strands of high-refractive index 
materials in air or in a homogeneous cladding, however, are far from optimal for detecting experimentally the effects 
of the novel term, because the maximum of the geometrical nonlinear coefficient is located, as a rule, far away in a 
region of strong normal dispersion, where bright solitons cannot exist. In fact, as has been anticipated in Ref. |3|, the 
new nonlinearity is visible only in presence of solitons, and its effects being nearly invisible in the normal dispersion 
regime. 

However, one is by no means restricted to the use of a homogeneous cladding around the high- index core. Here we 
explore the possibility of introducing the PhN into a silica-based photonic crystal fiber (PCF, [5]) with a triangular 
arrangement of the holes [see Fig. Ufa)]. We call such a design a photonic crystal fiber nanowire (PCF-NW). Such a 
structure has the additional advantage that relatively long PhNs can be supported by the robust PCF cladding Q. 
In this Letter we theoretically demonstrate that, by means of careful choice of parameters in the design of PCF-NWs, 
one can move the maximum of the geometrical nonlinearity inside the region of anomalous dispersion. This will make 
it much easier to observe the new geometrical nonlinearity experimentally, in that it causes considerable suppression of 
the soliton self-frequency shift (SSFS, see Ref. [7]) from the very start of the propagation. Thus, just a judicious choice 
of the geometry of the PCF cladding around the nanowire, which affects the dispersive properties of the waveguide 
and the way the mode profiles change with wavelength, is sufficient for the emergence of this completely unexplored 
optical nonlinearity, a feat that would be extremely difficult if not impossible to achieve with other simpler geometries. 

The fiber geometry that we propose in this Letter is shown in Fig. []Ja). It is made of a silica PCF cladding 
with a triangular lattice of air holes, with pitch A = 1.4 fim and hole radius R = 0.56 /im. The central core, of 
radius R c = 0.5 /im, is made of a high refractive index tellurite glass (T2 composition taken from Ref. [§[), which 
ossesses an estimated nonlinear coefficient ri2 ~ 4 x 10 -19 m _1 W _1 , almost 15 times larger than that of fused silica 
Its Raman gain spectrum h(uu) is shown in Fig. QJb), and has a large peak centered around 20 THz. In this 
waveguide, light confinement is provided by total-internal-reflection at the core-cladding boundary. Air holes in the 
rings modify the dispersion in such a way as to match our requirements, as will be explained shortly. Typically, only 
one ring of holes is sufficient to obtain the dispersive features described below. Our choice of core material and PCF 
parameters has been dictated by two conditions, which cannot easily be simultaneously met. The first requirement is 
that there must be a relatively large refractive index contrast between core and cladding, so that the magnitude of the 
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FIG. 1: (Color online) (a) Schematics of the proposed PCF nanowire geometry. The central nanowire core is made of tellurite 
glass (77TeO2-10Na2O-10ZnO-3PbO composition, T2 glass from Ref. [8(|), while the cladding is made of silica with a triangular 
lattice of air holes. Parameters are: pitch A = 1.4 /xm, hole radius R = 0.56 /xm, core radius R c = 0.5 /mi. Inset, the mode 
profile ||e(r_L)|| of the waveguide for A = 1.55 /im is shown, (b) Profile of the Raman gain spectrum for the T2 composition. 
Solid blue line: experimental data taken from Ref. [8[. Dotted red line: fit used in our simulations. 



longitudinal component of the electric field becomes appreciable, and in turn that the mode profile changes strongly 
with frequency. The second requirement is that the holes should considerably modify the GVD, in such a way that 
the maximum of the geometrical nonlinear coefficient is located inside the region of anomalous dispersion. However, 
the larger the refractive index contrast, the more the field is localized in the core, which makes the dispersion of the 
waveguide progressively more and more like that of a single rod surrounded by homogeneous silica, which has been 
proved to be non-optimal in Ref. A trade-off between the above conditions must be found. We have examined 
the dispersion of many high-index glasses for the core material (chalcogenide, bismuth, germanium-doped silica), and 
even some nonlinear liquids such as carbon disulfide, and systematically explored hundreds of specific parameters 
for the PCF cladding, but the above tellurite-based design seems to be one of the best solutions for our purposes. 
Such multi-glass hybrid waveguides can be fabricated by using the pressure-cell approach [10]. Melted tellurite glass 
is pressed under large external gas pressure into the holes of silica PCFs. The technique relies on the fact that the 
tellurite glass has a significantly lower softening point than silica, leaving the geometry defined by the silica host 
unchanged during filling [ll|. The problem of surface-induced glass crystallization in tellurite compound glasses has 
been solved in recent experiments and will be published elsewhere. Altogether, the pressure-cell technique provides a 
potential way to realize the proposed structures. 

For each frequency cj, the linear fundamental mode of the waveguide of Fig. QJa) has a normalized electric field 
profile given by e w (rj_), where are the transverse coordinates. The inset in Fig. [TJa) shows a contour plot of the 
norm | |e| | = [e^+e^+e 2 ] 1 / 2 for the fiber parameters given in the caption, corresponding to our representative PCF-NW 
design that we shall use throughout the paper. Due to the invariance of the PCF cladding under the Cq v symmetry 
group, the waveguide does not exhibit any birefringence, and it has a fundamental mode that is degenerate in the two 
orthogonal polarization states. The crucial point of our formalism is that one can describe the full ^-variations of e 
through the Taylor series 



Mrx) = V^(r,)f^V 



where Auj = uj — ujo is the frequency detuning from an arbitrary reference frequency cjo, and the quantity = 
ojQd^e UJ (r±)/du^ is proportional to the j-th frequency derivative of the mode profile. From now on, letters 

jhpv will be used for derivative indices. Following Ref. [3], one can rigorously prove that the equation governing the 
nonlinear light propagation of one of the two polarization states of the fundamental mode of the PCF-NW is given by 

id z Q + D(id t )Q + l jhpv &{id t )4> hpv = 0. (2) 

jhpv 

In Eq. ([2j), Q(z, t) is the electric field envelope, D(idt) = P(^o + idt) — /?(^o) * s ^ ne dispersion operator that encodes all 
information on the fiber GVD around ujo [l2[, /3 is the cj-dependent propagation constant of the fundamental mode, 
and & (idt) = [1 + (i/uJo)dt] [(i/uJo)dt] J is an operator that naturally contains the dynamics of the shock term at the 
zero-th order of the Taylor expansion (j = 0). The convoluted nonlinear fields used in Eq. (j2j) are defined as: 



, hpv(z t) ^ [(idt) h Q] WjWgj K-idt) v Q*])} 
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where symbol (g) is used to denote a time-convolution: A B = ^4(t — t f )B(t f )dt f =B A. In Eq. (j3]), = 
(1 — + Qh(t) is the nonlinear response function of the core, which is made of a tellurite glass following the design 

of Fig. Ufa), and includes the instantaneous Kerr [proportional to the Dirac delta 5(t)] and the non- instantaneous 
Raman [proportional to h(t)\ responses exhibited by the core material, being the relative importance between the 
two. In this Letter, from the experimental Raman gain of tellurite glass (solid blue line in Fig. DJb), taken from Ref. 
[H) we can extract a fit of h(uS) [dotted red line in Fig. DJb)] that we thus use in the numerical simulations. The 
dependence of R(t) can be safely neglected, since most of the energy is located in the core material, which is also 
much more nonlinear than the surrounding silica cladding. Note that in expression (j3]), (j) 000 = Q f R(t — t')\Q(t')\ 2 dt' 
gives the conventional, zero-th order convolution that is used in the GNLSE pH ]. 
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FIG. 2: (Color online) Linear and nonlinear data for the PCF-NW design of Fig. [T(a). (a) Blue solid and red dashed lines 
indicate respectively 70 and 71 versus wavelength, (b) Blue dots indicate the GVD of the waveguide, with zeros-GVD points 
located at A ~ 0.84 /im and A ~ 2.4 /im. Red solid line indicates parameter r = 7 1000 /7 0000 5 that appears in Eq. ©. 



The last ingredient in Eq. (|2j) contains the generalized nonlinear coefficients ji hpv , defined as 



y^ Ms -±J dr , x m x{r±) l _ 1, (4) 

(3) 

where Xxxxx(r±) is the third-order susceptibility, which is a function of the transverse coordinates in the waveguide of 
Fig. Ufa). Definition (j4]) is a generalization of the nonlinear coefficient commonly used in fiber optics [l2|, and takes 
into account the full vector nature of the field profile as well as its frequency variations. Such variations are at the 
very core of the new geometrical nonlinear it ies described here, since the Taylor series of Eq. (pQ) implies the existence 
of an infinite number of additional nonlinear terms that depend on progressively higher-order time derivatives of the 
envelope. The quantities ji hpv satisfy general symmetry relations that drastically reduce the number of independent 
nonlinear coefficients for each order of the derivative jH, G3 : ^ hpv = j vhp i = ^ phv = j vph ^ ^o hhh — ^hjhh ^ 
ryjhhj _ ryhjjh _ j n and (j3]), each field derivative is associated with a factor co^ 1 , which ensures convergence. 

Thus the physically meaningful nonlinear coefficients can be defined as ^ hpv = ^ hpv (ujQto)~^^ h ^ p ^ v \ where to is 
the input pulse duration. 

From the expression in Eq. (j4]) one can define the zero-th order nonlinear coefficient of the waveguide 70 = 7 0000 = 
7 0000 , corresponding to the conventional definition used in nonlinear fiber optics [12j. To first order in the Taylor 
expansion in Eq. (|2j) one can define the coefficient 71 = 7 1000 = 7 0100 = 7 0010 = 7 0001 = 7 1000 / (uJot ) , associated 
with nonlinear convoluted fields (j) hpv that contain only one time derivative of the envelope. Fig. [2fa) shows plots of 
70 and 71 versus reference wavelength for the fiber design of Fig. QJa). The fiber GVD is shown in Fig. [2^b) with 
blue dots. It is clear from this figure the well-known fact that 70 decreases monotonically for longer wavelengths [l5| . 
However, it is interesting to note that the geometrical nonlinear coefficient 71 initially increases, but then reaches a 
maximum near the infrared zero-GVD point of the fiber (located at A ~ 2.4 /im), and then tends to zero for even 
longer wavelengths. The close vicinity of max (71) to the anomalous GVD of the fiber is an atypical feature, that 
we have found only in a few very specific designs, including the one presented in Fig. QJa). The 'normal' situation, 
which is also true for circular PhNs surrounded by homogeneous media (such as, for instance, tapered fibers), is that 
max (71) is located well within the region of normal GVD [3[. 

We now show that the onset of geometrical nonlinearities leads to a strong suppression of the SSFS. In order to do 
this, we compare two direct numerical simulations, the first one obtained by truncating the Taylor expansion of Eq. 
(|2]) to the zero-th order (1 nonlinear convolution, corresponding to the conventional GNLSE) [Fig. 13(a)], the second 
one obtained by truncating the same sum to the second order (15 nonlinear convolutions in total), which thus takes 
into account the dominant terms of the geometrical nonlinearities [Fig. [3(b)]. It is easily seen by comparing Fig. [3(a) 
with Fig. [3(b) that the net effect of the additional nonlinearities is to reduce considerably the SSFS in the fiber. Thus 
a unique interplay between the geometrical nonlinearities and the SSFS takes place in properly designed PCF-NWs. 
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FIG. 3: (Color online) (a) Evolution of a to = 150 fs pulse in the waveguide of Fig. [TJa) according to Eq. ([2} truncated at 
the 0-th order, i.e. by using the conventional GNLSE. Second order dispersion length is Ld2 — 6 cm. (b) Same as (a) but 
truncating Eq. (|2} at the 2-nd order. The suppression of SSFS in the latter case is evident. Vertical white lines indicate the 
two zero-GVD wavlengths. (c,d,e) z-evolution of pulse duration r [to = r(0) = 30 fs], chirp C and soliton frequency shift f2 
according to Eqs. J3][7]), for aiV = 7 soliton, pump wavelength A = 1.4 /xm. Blue solid (red dashed) line refer to the case r = 
(r = 0.2) in Eqs. (loTfl) . 



It is possible to qualitatively understand the reason of the above SSFS suppression mechanism by using the so-called 
moment method p~2|, [H|. One assumes that after formation, each solitonic pulse does not change its functional shape, 
given by Q(t) = [Pt /r] 1/2 sech ([t - T]/r) e -^(t-T)-iC(t-T) 2 /(2t 2 ) ^ where T ig the tempora i delay of the solitonic 
pulse, Q is its frequency detuning from the reference frequency c^o, r is the soliton pulse width, C is the pulse 
chirp, fa = (d 2 /3/dcj) UJ = U j < is the second-order (anomalous) dispersion coefficient, P = N 2 Pq is the peak power, 
P = |/^2 1 (^oTo) - 1 is the fundamental (N = 1) soliton power, N is the soliton order and Tr = f_™ R(t)dt is the first 
moment of the Raman response. One can prove that the z-evolution of r, C and Q is then given by the following 
closed set of equations: 



dr 

dz 
dC 
dz 

dn 

dz 



T 

m\N 2 r 



1 



7T 2 toT 

8T R N 2 \(3 2 \ 1 

15t T S 



+ — (l + 4r) 
w 



(5) 
(6) 

(7) 



with the initial conditions r(0) = to, C(0) = Q(0) = 0. Higher-order terms in the dispersion and small terms 
proportional to C 2 have also been neglected for sake of clarity. Geometrical nonlinearities are parameterized by the 
coefficient r = ^ 1000 /^ 0000 5 the only one that appears in Eq. (j6]), shown in Fig. Efb). Note that the condition r > 1/4 
is necessary for the new nonlinearity to dominate the shock term in Eq. ((6]). Equations (|5]|7j) are written under the 
simplifying assumption that the GVD does not change during the soliton evolution, and that the pulse duration is 
longer than 100 fs, so that one can use a well-known approximate expression for the Raman term [12|- The term on 
the right-hand side of Eq. ©, proportional to Tr, is responsible for the constant SSFS along the fiber 0, The 
rate of this shift is very sensitive to the actual pulse width, since it is determined by r -3 , and it is always directed 
towards negative detunings, i.e. towards the red part of the spectrum. However, due to the action of the right-hand 
side of Eq. ©, the soliton acquires a small chirp even if its initial chirp vanishes. The slope of this chirp is initially 
negative, due to the initial condition r(0) = to, which makes the term proportional to Q/cjq dominant. Thus C < 
in the initial stage of propagation, which in turn leads to a pulse broadening due to Eq. (|5j), for which foC > 0. 
Finally, such broadening leads to a sharp suppression of the SSFS given by Eq. ([7j), due to the r -3 dependence of 
its right-hand side. This mechanism qualitatively explains the overall suppression of the SSFS due to the geometrical 
nonlinearity on soliton propagation observed in the direct numerical simulations of Eq. (j2j) . The long term behavior 
of the propagation can only be understood by numerically solving Eqs. (07j). I n Figs. E^c-e) we show the z-evolution 
of the quantities r, C and ft for parameters given in the caption. C undergoes several deep oscillations, and its sign 
is mostly negative throughout the whole propagation [Fig. [3fd)]. The pulse duration r undergoes similar oscillations 
(with a smaller magnitude), but overall it constantly grows [Fig. Efc)]. ^, however, is not too sensitive to such 
oscillations, due to the fact that its derivative never changes sign in Eq. (0) [Fig. [31(e)] . 

In conclusion, we have shown the emergence of a new type of nonlinearity in tellurite PhN embedded in a PCF 
cladding, which strongly depends on the geometrical parameters of the PCF design and the specific dispersion of 
the core material. Apart from its value as a new fundamental entity in nonlinear fiber optics, the existence of such 
nonlinearity shows that there is still much unexpected physics to unveil in complex PCF geometries. 
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